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ABSTRACT 


The MSC/NASTRAN finite element computer program and a Cray XMP 
computer were used to study the modal characteristics of a violin with the 
Stradivari shape. The violin geometry was modeled using an arcs of circles 
scheme with PATRAN, a finite element graphics pre/postprocessor program. 
The violin was modeled in-vacu and with free boundary conditions. Belly, 
back, sound post, bassbar, neck, bridge, tail-piece, strings, rib linings, end 
and comer blocks are the components of the model. Mode shapes and 
frequencies were calculated for free top and back plates, the violin box, and 
the complete violin system. Including the strings. The results obtained 
from the finite element technique were compared to experimental data from 
real violins collected by other investigators. 
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I. INTRODUCTION 


A. PAST AND CURRENT USE OF THE FINITE ELEMENT METHOD IN THE 

ANALYSIS OF STRINGED INSTRUMENTS 

In recent years there has been a vast increase in the application of the 
finite element method to diverse fields of the physical sciences. This 
explosion first reached the field of stringed instrument acoustics in 1974, 
with Schwab's finite element description of the modal shapes of a guitar 
top plate [Ref. 1], Richardson and Roberts [Ref. 2] continued the work on 
guitar plates in 1983, and Jenner [Ref. 3] is currently using the finite 
element method in his guitar research. 

In 1984, Richarson and Roberts [Ref. 4] and Rubin and Farrar [Ref. 5] 
applied the finite element method to the study of violin plate vibrations. 

The first descriptions of the mode shapes and frequencies of the violin 
box using finite elements were obtained in late 1984 by Richarson and 
Roberts [Ref. 6] and Knott [Ref. 7]. This study is simply a continuation of 
the previous research and has as its goal the modal description of the 
entire violin system. 
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B. THE PLATES STUDIED AS ENTITIES 


Before attempting to describe the modal characteristics of the finite 
element violin as a complete system the constituent parts must be 
studied as entities; then the physical knowledge of each part can aid In the 
attempt to explain the physical characteristics of the system as a whole. 
This study Implements the above technique. First, the top plate (see 
Figure I.O) and back plate (see Figure 1.1) of the violin were modeled with 
the finite element graphics pre/postprocessor language PATRAN . The 
natural frequencies for each plate were then obtained using 
MSC/NASTRAN*, a finite element analysis program (All analysis was 
performed on a CRAY-XMP* computer). Then the baseline dimensions and 
material properties of the plates were changed to discover how such 
variations affect the mode shapes and their frequencies, host of the 
effects of the variations were quite Intuitive and are not reported In this 
study. For example, Increasing the thickness of a plate shifts the natural 
frequencies up. Or, for another example. Increasing the across-grafn 
stiffness raises the frequency of those modes which have predominant 

*PATRAN Is a trademark of PDA Engineering, Santa Ana, CA, 92705. 

MSC/NASTRAN is a trademark of the MacNeal-Schwendler 

Corporation, Los Angeles, CA, 90041. 

CRAY-XMP Is a trademark of Cray Research Inc., Minneapolis, MN. 
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Figure 1.1. Back plate detail 



























































































bending action across the grain. Rodgers [Ref. 8] Is currently using the 
finite element method in a sensitivity analysis of the free violin plates to 
give a detailed description of how the mode shapes and frequencies are 
affected when various material properties In the plates are varied. 

C. THE VIOLIN BOX 

1. Modeling the Box 

The second stage of the study puts the belly (the top plate) and back 
together, Joined with the ribs, a luthier’s term for the side walls of the 
violin. The decision was made to include eight other components of the 
violin in this second stage: the soundpost; a dowel of wood attached 
vertically between the belly and back, the bassbar; a strip of wood 
attached to the underside of the belly, the two endblocks; blocks of wood 
attached to the Inside ends of the violin box for the purpose of structural 
rigidity, the four corner blocks; again solid blocks of wood attached to the 
'points’, or comers of the violin for the same purpose as the end blocks, 
and the rib linings; wood firring strips which strengthen the bond between 
the ribs and the belly and back (see Figures 1.2 - 1.5). 

2. The Analusls 

With the ’box’ modeled the second stage analysis could begin. The 

flrst.nine natural frequencies and mode shapes of the box model were 
12 
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Figure 1.2. Violin corpus 
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Figure I A Violin corpus - cutaway view 







Figure 1.5. Violin corpus - cutaway view 




calculated using the MSC/NASTRAN software. The mode shapes and 
corresponding frequencies were then compared to those obtained from an 
Impact hammer modal test conducted by Marshall [Ref. 9] of a real violin 
(neck, bridge and strings attached). An experimental modal test of a violin 
box sans attachments would have been desirable but It was found that 
Marshall's results confirmed the results of the finite element run In most 
respects. The mode shapes experimentally obtained were quite similar to 
those generated by the finite element run for those modes Involving the 
violin box alone. Of course the neck, bridge and string modes were absent. 
In addition the violin finite element model was in vacu so none of the air 
modes or their Influence on the wood modes was observed. The modal 
frequencies were also verified to be In the same general order and close to 
those obtained by Marshall. 

D. THE WHOLE VIOLIN 

The final stage of the study began with the modeling of the 
remaining components of the violin. With neck, bridge, strings and tall 
piece models added to the data base a complete violin system was realized 
(see Figures 1.6 - 1.8). 
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2. Analysis of the Whole Violin 


With the whole violin modeled the final step in the study was to run 
a normal mode analysis and once again compare the results to those 
obtained by Marshall's modal test on a real violin. This final analysis 
proved to be the most difficult portion of the study. All of the previous 
runs hod used conned solution sequences included in the MSC/NASTRAN 
software. Linear static and dynamic analyses of relatively simple 
problems can usually be solved using these preprogramed sequences of 
matrix manipulations. The addition of the strings and the requirement 
that they be under tension complicated the system to the extent that the 
rigid format solution sequences were no longer applicable to the structure. 
Specifically, off-axis deformations of the strings and large relative 
displacements between the string elements and those in the rest of the 
violin model introduced geometric non-linearities which must be included 
in the finite element model. 

The NASTRAN finite element analysis package is written entirely in 
DMAP (Direct Matrix Abstraction Programming), a high-level language for 
defining matrix operation algorithms. The geometric non-linearities in 
the violin system under tension were incorporated by using the DMAP 
language to alter the rigid format solution sequences. A brief discussion 
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of the DMAP language and the details of the alterations are discussed In 
the analysis portion of this paper. 

Once the solution sequence was altered the normal modes of the 
entire system were computed to 1,200 Hz. Thirty-nine mode shapes and 
frequencies were obtained and compared to those found by Marshall on a 
real violin. This time, as was expected, the agreement between 
experiment and the finite element simulation was even better than in the 
case of the box alone. 

E. SOURCES FOR THE INPUT DATA AND EXPERIMENTAL VERIFICATION 

The correct implementation of any numeric method requires that the 

input data be accurate, the methods be sound, and that the results be 

experimentally verifiable. As a result this work, like so much research, 

leans heavily on the efforts of others. Specifically, Eckwall [Ref. 10] and 

Sacconl [Ref 11] provided the geometry for the violin model. Eckwall's 

arcs of circles scheme for describing a Stradavarl violin saved hundreds of 

hours in data preparation. Haines [Ref. 12] and the United States 

Department of Agriculture [Ref. 13] was the source for most of the 

material properties. Marshall's recent study [Ref. 9], a modal analysis of a 

real violin, was an extremely valuable source for experimental 

verification and will be referred to often in this paper. Hutchen's work on 
22 



plate tuning [Refs. 14 and 15] provided the experimental basis for 
comparison of the finite element plate mode shapes with the mode shapes 
of real violin plates. 
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II. THE FINITE ELEMENT MODEL OF THE VIOLIN 


Ekwall’s work describes "II Cremonese 1715” ex Joachim entirely in arcs 
of circles [Ref. 10]. These arc lines are ideally suited as inputs to the 
PATRAN preprocessing graphics software. Once the outlines and archings of 
the violin box were described with the arcs, adjacent arcs were used to 
describe surfaces. The components of the box were then attached to the 
surfaces, resulting in a highly accurate numerical representation of the 
violin geometry. The finite elements were defined from the surface 
geometry and consist of triangular and quadrilateral shell elements for the 
belly, back, and tail piece, beam elements for the strings, and solid 
elements for the rest of the components. 

The material properties used for the analysis runs are listed in Table I. 
The data listed in Table I can be considered as a good starting point since it 
represents "average" values obtained from the wood handbook [Ref. 13] and 
from Haines' study [Ref. 12]. The belly, soundpost, rib linings, bridge, end 
and corner blocks are modeled as Sitka spruce; the back, neck, and tailpiece 
are modeled as Norwegian maple; and the bassbar is modeled as Douglas fir. 
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TABLE I 

MATERIAL PROPERTIES 

NECK, BACK, BELLY, RIBLININGS, 


TAILPIECE, RIBS 
(MAPLE) 

POST, BRIDGE,BLOCKS 
(SPRUCE) 

BASSBAR 

(FIR) 

El (Young'sModulus, 
PaX 109) 

17.000 

13.000 

12.400 

e 2 

1.564 

0.890 

0.621 

e 3 

0.731 

0.649 

0.621 

Gj (SheerModulus, 

1.275 

1.015 

0.794 

PaX I0 9 ) 




CM 

CD 

.1.173 

0.715 

0.868 

G 3 

0.187 

0.416 

0.967 

Kj 2 (Poisson's Ratio) 0.318 

0.375 

0.292 

•'13 

0.392 

0.436 

0.449 

*23 

0.703 

0.468 

0.390 

*32 

0.329 

0.248 

0.287 

*!21 

0.030 

0.034 

0.020 

*1 

0.019 

0.022 

0.022 

d (Density, kg/m 3 ^ 

800.000 

460.000 

506.00 
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The material orientations for the whole model are defined in Figure 2.0 
and are illustrated to avoid confusion when refering to the material 
properties Table. 


3/T 



R - Radial direction of wood grain. 

T - Tangential direction of wood grain. 

L - Longitudinal direction of wood grain. 

Figure 2.0. Material orientations 


The thickness of the belly was defined in accordance with Sacconi 
[Ref. 11 ] to be 3 mm except at the edges and area around the two ff holes 
(these holes can be seen in Figure 1.0), where it was 3.8 mm. In the violin 
model the top plate was modeled with a uniform thickness of 3.0 mm. The 
back was defined to the thickness of 4.5mm at the center, descending to 3.8 
mm on the first circles, to 3.5 mm on the second, to 2.5 mm in the upper 
lung and to 2.6 mm in the lower one. Figure 2.1 shows a countour map of the 
back thicknesses. 
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The neck dimensions were unavailable for the Stradivari violin, so the 
modeling of the neck was accomplished using neck dimensions from the 
author's violin. Hexahedral and wedge elements were used for the neck 
model. Neck details are Illustrated In Figure 2.2. 

For the tail-piece quadrilateral and triangular shell elements were used. 
The geometry of the tail-piece was obtained from the author's violin. 

The bridge model was created using the classical bridge form found In 
most violins. The model realization Is, of neccesslty, a simplification of 
the true geometry as can be seen In Figure 2.3. However, the obvious design 
objectives are retained. The thickness of the bridge varies from 5.0 mm at 
the feet to 0.5 mm at the top. 

The strings were modeled as beam elements. Beam elements were 
chosen instead of rod elements to allow the Inclusion of bending stiffness 
Into the string model, which the rod element does not consider. The 'G', 'D', 
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and 'A' string's material properties were modeled as gut cores with steel 
windings. The windings were modeled as a non-structural mass since their 
contribution to bending stiffness is assumed to be negligible. Solid steel 
was used for the 'E' string. Table II lists the dimensions and material 
properties for the strings. The number of beam elements for each string 
was determined from the desired frequency response of about 8,000 Hz. 


TABLE II 

STRING DIMENSIONS AND MATERIAL PROPERTIES 
STRING 


'G' 


CORE MASS(g) 0.110 

N.S.M.(g/m) 3.314 

TOTAL MASS(g) 1.200 

CORE DIAMETER(mm) 0.600 

E (Young's Modulus, Pa x 10 g ) 30.000 
v (Poisson's Ratio) 0.400 

d (Density, kg/m 3 ) 1200.000 

No. OF ELEMENTS 110 


'D' 

'A' 

'E' 

0.065 

0.034 

0.130 

1.352 

0.557 

0.000 

0.750 

0.260 

0.130 

0.460 

0.330 

0.270 

30.000 

30.000 

195.000 

0.400 

0.400 

0.280 

1200.000 

1200.000 

7700.000 

82 

46 

33 
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III. THE ANALYSIS 


A. REAL EIGENVALUE EXTRACTION TECHNIQUES 


The motion of the violin con be described by o set of equotions 
representing the bolonce of forces octing on the system. In the eigenvolue 
problem we are interested in the solution of the undamped and unforced 
system of equations: 

Mu + Ku = 0 (3.1) 


assuming a solution: 

u = ae 1wt (3.2) 

where a is a vector of displacement amplitudes, w is the natural frequency, 
and t is the time. Rearranging: 

(K -AH)a = 0 (3.3) 

where ?\= w 2 . 

The solution is nontrivial only if: 

det(K -?M) = 0 (3.4) 

A polynomial expansion of (3.4), where the values ofX which satisfy the 
equation are called the eigenvalues. 
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Each of the roots satisfies (3.4) so we have: 


(K - A 1 M)a 1 = 0 (3.5) 

where a^ is the vector of displacement amplitudes, the eigenvector, 
associated with the eigenvalue A j. An eigenvalue is obtained by solving 
(3.4) and then the eigenvector is determined by solving (3.5) with that value 
for A. 

In MSC/NASTRAN one of the available methods for extracting eigenvalues 
is the Givens method [Ref. 16]. For the Givens method the standard form of 
the eigenvalue problem is used: 

(J -Al)w = 0 (3.6) 


Where J is real and symmetric (as is the stiffness matrix K). Equation (3.6) 
is formed by decomposing the mass matrix, M, into Choleski factors C such 
that: 


M = CC T 

(3.7) 

Substituting (3.7) into (3.3) gives: 


(-ACC t + K)a = 0 

(3.8) 

So: w = C T a 

(3.9) 
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Substituting and rearranging gives: 

J = C _1 KC" t (3.10) 

and the standard form, equation (3.6) is realized. 

MSC/NASTRAN uses the Givens method to reduce the J matrix in (3.6) 
to tridiagonal form and then uses the Q-R method [Ref. 17] to obtain a 
diagonal matrix whose elements are the eigenvalues. The eigenvalues are 
then determined by the tridiagonal form of the eigenvalue problem. 

The computer runs used the above technique for eigenvalue and 
eigenvector extraction. Further explanation of the eigenvalue problem as it 
applies to MSC/NASTRAN can be found in the M5C/NA5TRAN Primer [Ref. 15]. 

B. TOP AND BACK PLATE ANALYSIS 

The eigenvalue extraction techniques covered in section A were used for 
the top and back plates. Appendix E, Figure E1, shows a simplified example 
of the data deck used. The initial runs on the top plate were without the 
bassbar and then with. Then the back plate was run. The analysis runs used 
free boundary conditions but rigid body modes are not included in the tables. 
Table III lists the natural frequencies, Figures A1 - A6 illustrate the mode 
shapes for the top plate with bassbar and Figures B1 - B6, the back plate. 
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TABLE III 

NATURAL FREQUENCIES FOR FREE TOP AND BACK PLATES 


Mode 

Top(Hz) 

Top/w/Bar (Hz) 

Back 

1 

67.0 

88.0 

97.0 

2 

115.0 

161.0 

170.4 

3 

200.0 

255.3 

268.2 

4 

202.3 

272.7 

268.3 

5 

274.0 

330.4 

353.0 

6 

- 

422.0 

413.7 

7 

- 

436.0 

453.7 

8 

- 

494.1 

558.0 

9 

- 

617.4 

632.4 
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C. THE VIOLIN BOX 


The violin box used the same executive routine directives as the free 
plates for its eigenvalue analysis runs. Of course the bulk data deck, which 
specifies the model geometry and material properties, was much larger in 
that it included both plates and the additional components making up the 
violin box model. Table IV lists the nine extracted natural frequencies and 
their descriptions. Figures Cl - C18 illustrate the corresponding mode 
shapes. Note that since the displacement amplitudes are mass normalized 
for each eigenvalue and are therefore meaningless for comparison of 
amplitudes among distinct modes, absolute displacements are not shown on 
any of the contour diagrams in this paper. 

D. THE WHOLE VIOLIN 

As mentioned in the introduction the transition from the violin box model 
to the model of the violin with all of its components was a difficult one. 

The greatest difficulty was in the fact that the finite element “strings" had 
to be streched a large amount relative to other strains in the violin. Large 
errors in strain energy were experienced when the complete violin system 
under string loading was programmed for a statics analysis. The rigid 
formats available in the NASTRAN software had to be modified, and the 
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TABLE IV 

NATURAL FREQUENCIES FOR F.E. CORPUS 


Mode Freq.(Hz) Description (IP = in phase, OP = out of phase) 


1 462.9 

2 478.9 

3 573.9 

4 682.5 

5 736.2 

6 749.6 

7 810.9 

8 899,1 

9 968.5 


Vertical translation of C bouts. Top and 
back IP. Left and right sides OP. 

Vertical bending of corpus. Top plate 
transverse, back plate longitudinal. 

Vertical bending of corpus. Top plate 
longitudinal, back plate transverse. 

Top and back plates OP. Essentially a 
lengthwise dipole. Back plate upper 
hump extends through the center area 
and resembles a ring mode. 

Dipole action in lower bout of top plate. 
Breathing in back plate, strongest in 
lower bout. 

Dipoles in lower bouts of both plates. 
Strong torsional action of the corpus. 

Similar to mode 4, a lengthwise dipole, 
except that top and back are IP. 

Double dipoles in top and back plates. 

Top and back OP. 

Top exhibits three “poles" in lower bout. 
Back plate shows ill-defined dipoles in 
upper and lower bouts. 
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DMAP language used to tailor the solution sequence to the problem. The 
following is a step by step explanation of the process: 

1. Assemble the stiffness matrix for the violin box with the neck 
included, but without the bridge/tail piece/strings system. A 
static analysis, SOL 61, was run with an alteration in line 355 of 
the preprogrammed solution sequence (see Fig. E2 for the data deck 
listing). The ALTER statement simply provides the capability 
to extract the stiffness matrix for the follow-on runs (named 
KGGW in Fig. E2) 

2. Perform static condensation of the above stiffness matrix 
"eliminating" all degrees of freedom except at those 14 locations 
where: the strings attach to the neck (4), where the bridge feet attach 
to the top plate (8), and where the two gut cords from the tail piece 
attach to the body (2). MSC/NASTRAN's substructuring capability 
was utilized for the condensation. The 5ESET cards in the data deck 
(see Fig. E2) define the substructure, in this case all of the nodes in 
the violin box/neck that are to be omitted from the KGGW matrix. In 
other words, the stiffness of the entire box/neck structure is defined 
for each unit stiffness of each of one of the 14 points, effectively 
“eliminating" all degrees of freedom except at those 14 points. 

3. Run an iterative, non-linear static solution sequence for the 
bridge/tail piece/strings system (see Fig. E3) using the 84 X 84 
KGGW matrix (14 nodes times 6 degrees of freedom per node) matrix 
obtained in step 2 as the boundary conditions. The stiffness matrix 
KGGW is essentially the stiffness of the neck and body that the string 
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ends and bridge feet would experience. Problems of large relative 
strains are avoided since the box/neck system is present only as a 
boundary condition (strains are not calculated for the box/neck). The 
first iteration uses the rigid format for a static analysis in which the 
strings are subjected to differing amounts of thermal stress. In 
essence, the strings were “cooled down" to shrink them and bring 
them to the required tensions. The desired tensions were those 
that would place the natural frequencies of the strings as close as 
possible to the tones found on a real violin, that is, G, D, A, and E. The 
following subcases (or iterations) of the run successively add the 
geometric non-linear terms and compute the stiffness matrix 
for each added term until an equilibrium of forces exists: 

F = P + Q (3.11) 

where F is the vector of element forces, P is the vector of applied 
thermal loads, and Q is the vector of forces of constraint. The 
string/bridge system stiffness matrix is now modified for the stress 
of the thermal loads on the strings and is saved for the next step. 

4. The final run puts the violin system back together (see Fig. E4). The 
string/bridge system, now in equilibrium, is added to the box/neck 
system. Then the rigid format solution sequence computes the 
normal modes of the entire violin model. 

Table V lists the natural frequencies obtained and a description of each 
mode. Figures D1 - D74 illustrate the mode shapes with contour diagrams 
and wire mesh deformed plots. Listings of the executive and case control 
decks for the analysis runs are included in Appendix E. 
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TABLE V 

VIOLIN SYSTEM NATURAL FREQUENCIES 


Mode 


number 

Frequency (Hz) 

Description (OP = out of phase, IP = in phase) 

1 

191 

Fingerboard bending and tailpiece 

deflection. 

2 

201 

Fingerboard bending and tailpiece 

deflection. 

3/4 

212/213 

G1 

5/6 

312/313 

D1 

7 

363 

Maximum motion in C bouts. Top and back 

OP. Left and right sides OP. 

8/9 

421/423 

G2 

10 

439 

First bending of corpus. Top plate 

transverse, back plate longitudinal. 

11 

458 

Same as mode 7 except top and back IP. 

12/13 

486/487 

Ai 

14 

494 

First fingerboard torsion. 
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TABLE V (CONT) 


Mode 

number Frequency (Hz) Description 


15 536 First bending of corpus. Top plate 

longitudinal, back plate transverse. 
Strong relative displacement on top plate 
in the thickest area of the bassbar 
inhibits the longitudinal bending 
somewhat. 


16 

568 

Fingerboard torsion. 

17 

618 

A combination of G3,D2, and El. 

16/19 

629/632 

D2 

20 

640 

G3 

21/22 

648/652 

El 

23 

656 

Another combination string mode. 

24 

683 

Relatively strong breathing action in 


upper bouts of both plates. First 
mode where bridge moves vertically 
rather than rocking back and forth. Top 
and back plates are OP. 
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TABLE V (CONT) 


Node 

number Frequency (Hz) Description 


25 

693 

Tailpiece twisting'. 

26 

727 

Dipole in lower bout of top plate. 

Relatively small breathing action of OP in 

lower bout of back plate. 

27 

771 

First torsional mode of corpus supported 

by a strong dipole in lower bout of back 

plate. The top plate lower bout shows a 

relatively low amplitude dipole. 

28 

812 

Strong IP breathing of lower bouts in top 

and back plates. Essentially a lengthwise 

dipole. 

29 

856 

G4 

30 

885 

Double dipoles. Top and back plates OP. 

31 

902 

G4, D3, and A2 combined string modes. 

32 

908 

Top plate exhibits three 'poles' in lower 


bout. Back plates shows dipoles in upper 
and lower bouts. 
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TABLE V (CONT) 


Mode 

number 


33 

34 

35 

36 

37 

38 

39 


Frequency (Hz) Description 


961 D3 

971 Similar to mode 32 in that it is of 

opposite phase except in the lower bout of 
the top plate. 

1000 A2 

1083 Three poles in lower bouts of both plates. 

Top and back OP. 

1092 Second bending mode of fingerboard/neck. 

Three poles in upper and lower bouts of 
both plates. Top and back IP. 

1121 Similar to mode 37 except no fingerboard 

action. 

1155 Back plate; three poles in lower bout, two 

in C bout, and two in upper bout. Top 
plate; similar except C bout area 
translates vertically. 
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IV. CQMBARI SQM JO. EXP.EB1MENT 


A. GENERAL COMMENTS 

The finite element method is a numerical technique that cannot stand 
alone; this is especially so when applied to a complex structure such as the 
violin. Ideally, as each subsystem is simulated using finite elements, a 
parallel verification of the model is accomplished using experimental 
methods. For the violin, an exhaustive verification would involve: 

1. Top and back plate. 

2. Bassbar. 

3. Top plate with bassbar. 

4. Violin box without sound post. 

5. Violin box with sound post. 

6. Neck. 

7. Violin box with neck. 

8. Bridge. 

9. Tail piece. 

10. Strings with rigid mounts. 

11. Whole violin. 

In the ideal case the material properties and precise geometry of each 
component used for experimental verification would be known and would 
constitute the data base for the finite element modeling process. 

In the present study steps 1,3, 5, and 11 were the only ones which were 
compared to experimental results. The experimental data was obtained 
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from Marshall [Ref. 9] and Hutchens [Refs. 14 and 15]. There exists no 
geometry or material property data for Marshall's violin SUS *295 or 
Hutchen s violin plates so the examination of the results in the present 
study must be done with the realization that natural frequencies and mode 
shapes were not expected to match exactly. A basic assumption to be tested 
in this study is that there exists a fundamental commonality of mode shapes 
among all but the most poorly constructed violins. This assumption has 
been shown to hold in the case of free violin plates, at least among the 
lower modes of vibration [Ref. 14]. That is, even plates with differing 
geometries and material properties show a clear similarity in the ordering 
and shape of their vibration patterns. That all violins sound like a violin 
(some sound better than others!), is the intuitive and rather obvious basis 
for the mode shape similarity assumption. 


B. TOP AND BACK PLATES 

The free top and back plate mode shapes of the finite element analysis 
were found to compare well with experiments done using Chladni patterns 
and holographic interferometry. Figure 4.0 is on example of holographic 
pictures of a back plate [Ref. 15] (see Figures B1 - B5 for the corresponding 
mode shapes obtained using the finite element method). 
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C. THE VIOLIN BOX 

Experimental data on a violin box was not available but the author 
thought It wise to attempt a comparison of the finite element box to that of 
a complete violin modal test. Then the transition to the complete finite 
element model would possibly be less confusing. Table VI compares the 
finite element box to the violin modal test results conducted by Marshall 
[Ref.9], SUS*295 refers to the particular violin used by Marshall. Bear in 
mind that the finite element model does not Include the neck, bridge, tail 
piece or strings, so the mode density of the computer model Is much less 
than that of the real violin. 
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TABLE VI 


MODE COMPARISON BETWEEN F.E. CORPUS AND SUS *295 


MODE * FREQ (Hz) 

F.E./SUS*295 F.E./SUS*295 

1/10 462.9/410.5 

2/12 478.9/466.1 

3/16 573.9/574.1 

4/19 682.5/656.1 


5/22 736.2/690.8 

6/24 749.6/760.3 

7/28 810.9/882.0 

8/29 899.1/930.1 

9/31 968.5/1087.0 


COMPARISON OF MODE SHAPES 


Similar. 

Similar. 

Similar. 

That the two modes are related is 
suggested by the presence of similar 
motions in the back plates. The top 
plates would not be expected to 
have the same motion since the one 
has the bridge as a major factor in 
the mode (it translates vertically) 
and the other has no bridge at all. 

No difference except breathing 
action is weaker in back plate 
of the finite element model. 


Similar. 

Similar. 

Similar except that in the 
finite element model the lower 
bouts are OP. 

Similar. 
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D. THE WHOLE VIOLIN 


Table VII details the comparison between the finite element model, 
complete with strings, etc., and 5US*295 used by Marshall for his modal 
test. The complete violin model compared reasonably well with the real 
violin in mode shapes, frequencies and relative order of the modes as can 
been seen by studying Table VII. The source of the discrepancies could be 
due to a number of factors, the most obvious being that the degree to which 
the material properties and geometries of the finite element violin and the 
real violin match is unknown. Any firm conclusions about the origins of the 
discrepancy between the finite element model and the real violin would 
require the actual dimensions and material properties of violin SUS-*295. 
The following are other causes of disagreement between mode shapes and 
frequencies between the finite element model and SUS -*295: 

- In the finite element model the strings were not tuned to the exact 
frequencies used on a real violin. For the 3 string the model was off 
by +16 Hz; the D string, +18 Hz; the A string, + 46 Hz; and the E 
string, - 7 Hz. The differences in tuning could have caused some of 
the discrepancies in mode shapes and ordering of the modes. Fixing 
the tuning problem will merely require an adjustment of the thermal 
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loads applied to the string model but, unfortunately, this was not 
possible prior to publication. 

The air modes of the real violin can be expected to couple to varying 
degrees with the structural modes of the violin. The model does not 
contain "air" elements. The air modes are absent from the finite 
element model and any effect they may have had on the structural 
modes is a matter for future research. 

A real violin is glued together in most places. In the finite element 
model there is no attempt to model the glue joints. 

Pressure from the plates and strings hold the soundpost and bridge to 
the real violin body. In the finite element model most of these forces 
exist (not the stress mounting of the soundpost in the violin box), but 
the components are "attached" to the plates and the bridge. In other 
words, there is no provision for slipping in the case of the strings or 
decoupling from the plates in the case of the soundpost. 
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TABLE VI! 

COMPARISON BETWEEN FINITE ELEMENT MODEL AND SUS *295 


Mode 

number 

F. E./5US*295 


1.2/1,7 


7,11/10 


10/12 

15/16 

16/11 

17/13 

24/19 


Frequency (Hz) Comparison 

F.E./SUS*295 


191,201 
/119,303 


363,458/410 


439/466 

536/574 

568/435 

618/471 

683/656 


In the F.E. model, the vertical 
deflection of the tailpiece found in 
mode 1 of SUS*295 and the first 
vertical bending of the neck fingerboard 
found in mode 7 of 5US*295, are both 
combined in mode 1. Modes 1 and 2 in 
the F.E. model are dual modes 
F.E. mode 7 is coupled with a sideways 
deflection of the tailpiece. F.E. mode 
11 matches 5115*295 mode 10 in phase. 
Similar. 

Similar. 

Similar. 

Similar. 

An interesting point here is that in the 
F.E. model the sound post is not stress 
mounted to the box. The difference in 
mode shapes between the F.E. model and 
SU5*295 seems to originate from this 
difference between the model and 

reality. 
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TABLE VII (CONT) 


26/22 

727/691 

Similar. 

27/24 

771/760 

Similar. 

28/28 

812/882 

Similar. 

30/29 

885/930 

Similar except that in the F.E. 



model the lower bouts are OP. 

32/31 

908/1187 

Similar. 


* For higher mode numbers the F.E. model and SUS*295 do not compare 
well. For example, all higher modes in the F.E. model have a larger number 
of poles. 
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V. CONCLUSION AND FUTURE WORK 


Since this paper is little more than a catalogue of a finite element 
violin's mode shapes more questions are raised than answered. The 
purpose of the study was to provide this preliminary catalogue and to 
lend support to the idea that the mode shapes of a violin, complete 
with the string system, can be described using the finite element 
method. It should be kept in mind that the method's real utility lies 
not in mere reproduction and description of physical reality but in 
sensitivity analyses and design optimization. However, before 
research can progress to violin design optimization still more will be 
required in the effort to accurately describe the violin. Further work 
will first concentrate on gaining a more complete data base of real 
violins and components, with known geometries and material 
properties. A coincident finite element analysis can then be expected 
to agree with experimental modal tests of the real violins and 
components. Different violins can be tested and modeled and the 
mode shapes compared, perhaps leading to an understanding of what 
differentiates a “good'' sounding violin from a "poor" one. 
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The incorporation of air elements into the model is a long-term but 
necessary goal for further research. The extent of the air mode's coupling 
effects with the structural modes could be easily determined through the 
study of simple finite element models. Radiation and string exitation 
studies using finite element models are also promising areas for further 
work. 

After developing the finite element model of the violin the author found 
the model to be nearly as full of mystery as a real violin! A more methodi¬ 
cal, step-by-step analysis should clear up many of the grey areas. The 
finite element method is an extremely powerful descriptor of physical 
reality and this author believes that much of the violin “puzzle" will be 
solved as research continues. 
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APPFNDIX A 


TOP PLATE MODES OF VIBRATION 


The contour figures are drawn to show motion as it 
would appear if viewed from above the plate, with plus 
signs representing motion up, or towards the viewer, and 
minus signs representing downward motion. 

The darker contour lines represent lines of near zero 
motion. 














Figure A3. Mode 3 Figure A4. Mode 4 
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Figure A5. Mode 5 Figure A6. Mode 6 




















AP P ENDIX 


BA C K PL ATE MODES QF V IBRATION 


The contour figures are drawn to show motion as it 
would appear if viewed from above the plate, with plus 
signs representing motion up, or towards the viewer, and 
minus signs representing downward motion. 

The darker contour lines represent lines of near zero 
motion. 


57 




58 






59 


Figure B3. Mode 3 Figure B4. Mode 4 
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A PPEND IX. C 


VIOLIN CORPU S HO PES OF V I BRAT I O N 


The contour figures are drawn to show motion as it 
would appear if viewed from above the plate, with plus 
signs representing motion up, or towards the viewer, and 
minus signs representing downward motion. The back 
plate is seen as if the top plate was transparent. 

The darker contour lines represent lines of near zero 
motion. 
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Figure C7. Mode 4 - top 682.5 Hz Figure C8. Mode 4 - back 
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Figure Cl 1. Mode 6 - top 749.6 Hz Figure Ci2. Mode 6 - back 

























APPENDIX D 


COMPLETE VIOLIN MODES OF VIBRATION 


Most of the figures in this appendix are in groups of four 
figures per mode. The first plot is of the deformed mesh 
as it would appear to the viewer looking at the top, or 
belly of the violin. The second plot is a contour diagram 
of the top plate viewed from above the instrument. The 
third plot is of the deformed mesh as it would appear to 
the viewer looking towards the back of the instrument, 
from below the violin. The fourth plot is a contour 
diagram of the back plate viewed from above the 
instrument as if the top plate was transparent. 

For the contour figures plus signs represent motion up, 
or towards the viewer, and minus signs represent 
downward motion. 

The neck, bridge and strings are not shown in the contour 
figures for purposes of clarity. 

The darker contour lines represent lines of near zero 
motion. 

The strings are not shown in the deformed mesh plots for 
purposes of clarity. 

For the string mode deformed plots, the rest of the finite 
element model is not shown. 
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536 Hz 

Figure D19. Mode 15 



































683 HZ 

Figure D27. Mode 24 







































Figure 042- Mode 28 



Figure D43. Mode 28 - back 
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1083 HZ 


Figure D59. Mode 36 



Figure D60. Mode 36 - top 
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APPENDIX E 


NASTRAN PROGRAM LISTINGS 


The following program listings are of the NASTRAN case 
control, executive routine, and partial bulk data decks 
for the violin model. 
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For the top and back plates a rigid format solution sequence was used. 
Figure El is a partial listing of the executive, case control, and bulk data 
decks. Comments are shown in light face type. 


SOL 3 Rigid format for normal modes analysis. 

CEND Designates the end of the executive deck. 

METHOD - 1 Set ID directive. Refers to El SR, 1 below. 

SET 1 = ALL All grid points (nodes) are in set one. 

BEGIN BULK 

EIGR, 1, MGIV, O.O, 200.0,22,,, I.-8,+EI 
+El,MASS Set ID 1 real eigenvalue extraction using 

the modified Givens method. Roots from 
0 to 200 Hz with eigenvectors mass 
normalized. 22 roots estimated. 

GRID, 1„0.,0.,0. Node number and x, y, z coordinates. 

• 

GRID, 1000„20.,30., 15. 

CTRI A3,1,1,6,13,7 Element type. 


CTRI A3,1900,1,258,259,260 

PSHELL, 1,1,0.3,1 „ 1 Property card. 

MAT8,1,13.* 10,8. + 9,.375,3.+9,9.-*9,.7+9,.46 

ENDDATA 


Material 
property card. 


Figure El. Normal modes analysis example data deck 


It should be emphasized that Figure El is apartial listing of the 
executive and case control decks. Refer to the Finite Element Primer 
[Ref. 15] for the detailed information needed to assemble the data deck. 




Figure E2. Equivalent Stiffness Matrix Data Deck 
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33-1) ISOTROPIC I'l'.PLE-LIJCE MATERIAL FOR NECK TETRA ELEMENTS 
MAT 1 ,6,17.5+9. , .310,0.0-4,1). ,0. 

PS0L11), 15,6 



OSTEEL FOR TAIL PIECE 

MAT1,4,19.5+10,,0.20,1.B6-3 

PS11ELL, 14,9.4. ,9, ,9 

SCUT CORE FOR G, D, AND A STR1NCSIPREAMI6-C, 17-D, 10-A, 20-TAIL) 
MAT1,0,3+10 ... 4,1.2-3,6.-6,0. 

PBEAi!, 1 9,7,57.3-3 .261 .3-6.261.3-6, , 522.6-6 
PDF.AM , I 6,0,202.0-3,6.320-3,6.323-3, , >2.636-3,2.692-3 
PDF.AH,17,0, 166.0-3,2.193-3.2.1 93-3, ,9.306-3,1.352-3 
PlSEAll, 111,0,25.5-3,501 .7-6,501 .7-6, , 1 163.-6,0.337-3 
PBEAI'l .20.0,3.2,434.3-3,404.5-3 ,.969.-3 

srnt ron bassd/.r 

M.AT9,3,12.697h <> ,400.65+6,396.63+6, , , ,963.59+6 ,+MlOO 
+.1190,205.5„+6 , , . ,71 1 .07+6 , , , , +H200 



Figure E2(Cont.) Equivalent Stlfness Matrix Data Deck 
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Figure E3. String Analysis Data Deck 
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Figure E3(Cont.) String Analysis Data Deck 
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SOL 63 
TIME 120. 


SUBTITLE = CREIlOriESE 1715 

LABEL = RESTART VIT1I STRING DIFFERENTIAL STIFFNESS 

ECHO = NONE 

SEALL = ALL 

SPC = 1 

DYNRED = 1 

NCTilOD = 1 

BEGIN BULL 

PAIIAM, MAXRAT 10,1. +20 

SPC1,1,0,99001.THRU.99150 

SPOINT,99001.THRU,99150 

ASET1,0,99001,THRU,9 V150 

QSETI.0.99901,TURU,99150 

DYNRED,1.1200. 

E1GR,1,GIV,0.,1200.,,,,1.-8,+El 
+E1,NAX 
S 

*READ DMIG 


Figure E4. Restart With Differential Stiffness 
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